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In this reply we show that the criticisms raised by J. Noronha are based on a misapplication of the 
model we have proposed in [A. Jaouadi, M. Telmini, E. Charron, Phys. Rev. A 83, 023616 (2011)]. 

Here we explicitly discuss the range of validity of the approximations underlying our analytical 
model. We also show that the discrepancies pointed out for very small atom numbers and for very 
anisotropic traps are not surprising since these conditions exceed the range of validity of the model. 


The study of the critical temperature T c associated 
with the condensation of dilute Bose gases has attracted a 
lot of attention in the last two decades [1-3] . In this con¬ 
text, our work published in Physical Review A in 2011 [4] 
is the continuation of a previous study devoted to Bose- 
Einstein condensation (BEC) in blue-detuned Laguerre- 
Gauss optical traps [5]. Two years after these theoretical 
investigations, BEC was finally observed experimentally 
using a hollow laser beam [6] , confirming our prediction 
that condensation could be reached in such an experi¬ 
mental configuration. 

In this series of two papers, Ref. [4] specifically focused 
on the influence of the external trapping potential on 
T c beyond the thermodynamic limit, for a finite parti¬ 
cle number N. It is well known that the value of T c 
is influenced by the different nature of the condensation 
process in homogeneous vs. inhomogeneous gases [7, 8] . 
In harmonic traps, Bose-Einstein condensation is a local 
phenomenon in the sense that, initially, it only affects 
the atoms located near the trap center. On the contrary, 
in the case of a homogeneous gas, the phase transition 
is dominated by long-range atomic correlations. Investi¬ 
gating the crossover between these two limits of homo¬ 
geneous and inhomogeneous potentials was one of the 
motivations for our study and we used power-law traps 
to bridge this gap. 

With this aim, we have used the local density approx¬ 
imation (LDA), which assumes that the gas can be con¬ 
sidered locally as homogenous, with a local chemical po¬ 
tential expressed as 

Tiocaiip) = /I- V(r), (1) 

where V (r) is the external trapping potential and /i is the 
chemical potential at the center of the trap. In typical ex¬ 
perimental conditions using harmonic traps, it is known 
that this approximation yields very accurate results [1]. 
As mentioned in our original article, it is also known 
that this approach fails in the fully homogeneous limit 
corresponding to cold atoms trapped in a flat box [7]. 
For this limit, obtaining quantitatively accurate descrip¬ 
tions requires more elaborate theoretical methods. Our 


study based on LDA is nevertheless interesting for two 
reasons. First, it is an approach from which analytical 
results can be derived. It therefore provides a useful first 
approximation, known to be valid in usual experimental 
configurations with harmonic traps, against which one 
can compare more sophisticated theories. Second, due 
to its simplicity, it leads to a better understanding of the 
crossover between the two limits mentioned previously, of 
homogeneous and inhomogeneous potentials. This type 
of insight is more difficult to gain from other theories, 
which do not necessarily provide analytical results. 

In our calculation of T c , the thermodynamic sums over 
discrete quantum states are replaced, using the LDA, by 
integrals in phase space. This approximation, which is 
very common [1], is of a semi-classical nature. To be 
valid, it requires the thermal energy k B T c to far exceed 
typical quantum state energy differences [3] . In the case 
of an isotropic harmonic trap, this is simply written as 

k B T c » hco , (2) 


where oo denotes the trap angular frequency. In addition, 
in this specific case the critical temperature for Bose- 
Einstein condensation of an ideal gas in the thermody¬ 
namic limit is given by 


T o = to ( N V 
C k B VC(3 )J 


( 3 ) 


where N denotes the total number of particles and ((s) 
is the Rieman zeta function. 

The approach derived in [4] is an approximation and it 
is therefore important to evaluate its domain of applica¬ 
bility. In the case of an isotropic harmonic trap this task 
is relatively trivial. It requires to make an arbitrary but 
reasonable choice for the range of validity of Eq. (2) . In 
the following, we consider that the discrete nature of the 
trapped states can be considered as a continuum when 


k B T c > 20 Hco. (4) 

Using Eqs. (3) and (4), we obtain the following limit of 
validity for our approach 

N > 20 3 C(3) ~ 9600. 


( 5 ) 
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FIG. 1. (Color online) Rescaled temperature T c /T° as a func¬ 
tion of the logarithm of the particle number. The dashed 
blue line shows the critical temperature calculated at first or¬ 
der and the solid red line shows the result of our generalized 
LDA formula [4], which includes higher-order corrections. The 
temperatures T 0 .i%, To. 5 % and T\% at which the condensate 
fraction reaches 0.1%, 0.5% and 1% are shown in the three 
dotted green lines (see the three labels in the graph). 



FIG. 2. (Color online) Condensate fraction as a function of 
the rescaled temperature T/T°. The curve seen on the left 
with solid circles is for N = 10 4 and the curve on the right 
with solid squares is for N = 10 5 . The condensation temper¬ 
atures calculated at first order and with our generalized for¬ 
mula are shown as dashed blue and solid red lines respectively. 
The green area represents a typical range where the conden¬ 
sation temperature can be detected experimentally [9, 10]. 


This simple estimation shows clearly that the LDA can¬ 
not be used with very small particle numbers. This lim¬ 
itation is known since a long time [7] . In an isotropic 
harmonic trap, one can therefore expect LDA to hold 
for N ^ 10 4 . From this simple evaluation, one can al¬ 
ready conclude that the results presented by J. Noronha 
in Fig. 1 of his comment for values of N as low as N = 100 
are very far outside the limit of validity of our model, 
and they are therefore necessarily misleading. One can¬ 
not draw any serious conclusion from a comparison per¬ 
formed in this limit of very small particle numbers. 

Now, in the case of non-harmonic traps, the limit of 
validity of our approximation has to be re-evaluated. In 
a disk-shape or in a cigar-shape harmonic trap, the in¬ 
equality (2) becomes 

ksT c ^> sHuj, ( 6 ) 


where s denotes the anisotropy parameter defined in the 
comment of J. Noronha. In this case, the condensation 
temperature in the thermodynamic limit is obtained as 


1? = s r ‘ 


hu> / N \ 3 

^ VC(3 )) 


( 7 ) 


where n = 1/3 for a disk and n = 2/3 for a cigar. The 
validity criterion (5) then becomes 

N > (20s (1 " n) ) 3 C(3), (8) 

or equivalently 


s < 


1 

20 


O’l 


( 9 ) 


We now apply Eq. (9) to the two cases studied by 
Noronha in Fig. 2 of his comment, with N = 10 5 . The 
range of validity of our approximation, as defined by 
Eq. (9), is s < 3.2 for a disk and s < 10.4 for a cigar. 
We can clearly see in Fig. 2 of this comment that the 
estimations based on our simple analytical model start 
to deviate from an exact numerical calculation (denoted 
as T 0 i% in this comment) when s > 3 in the case of a 
disk (upper plot) and when s > 10 in the case of a cigar 
(lower plot). This is in good agreement with the expected 
range of validity of our approximation. Clearly, one can¬ 
not draw any serious conclusion from comparisons per¬ 
formed in the limit of very anisotropic traps. Conversely, 
in the range where our approximation is justified, that 
is for a relatively small anisotropy, we can conclude from 
the calculations performed by J. Noronha that the results 
given by our simple analytical formula agree rather well 
with other more involved numerical estimations of T c . 

Let us now look at the results presented by Noronha 
in more details. In the case of an isotropic trap, we re¬ 
produce in Fig. 1 the variation of the rescaled conden¬ 
sation temperature T c /T® as a function of log 10 (AT) in 
the range of validity of our model, i.e. for N ^ 10 4 . 
The dashed blue line shows the condensation tempera¬ 
ture calculated with the well known first-order finite size 
correction, while the solid red line shows the result of our 
generalized LDA formula [4] , which includes higher-order 
corrections. The temperatures T 01 %, T 05 % and T x % at 
which the condensate fraction reaches 0.1%, 0.5% and 1% 
are shown as different dotted green lines. These three ad¬ 
ditional temperatures have been evaluated directly from 
thermodynamic sums over the discrete energy levels of 
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the trap. We see in Fig. 1 that the first order formula 
lies between Ti% and T 0 . 1 %, while our result lies between 
T 0 5 y 0 and T 0 !%. From this graph, one can only conclude 
that the two analytical approximations give similar and 
realistic results, which converge to the same thermody¬ 
namic limit for large N values. 

In addition, in Fig. 2 we show the condensate fraction 
as a function of the rescaled temperature T/T°. The 
curve seen on the left with solid circles is for a num¬ 
ber of atoms N = 10 1 2 3 4 and the curve on the right with 
solid squares is for N = 10 5 . The condensation temper¬ 
atures calculated with the first order correction formula 
and with our generalized LDA formula [4] are shown as 
vertical dashed blue and solid red lines respectively, us¬ 
ing the same color code as in Fig. 1. We see here, again, 
that our generalized formula yields a slightly higher con¬ 
densation temperature than the one obtained with the 
first order approximation. It is therefore associated with 
a smaller condensate fraction. 

In typical experiments, the condensation temperature 
is measured by fitting the bi-modal distribution of the 
gas. Using such a scheme, condensate fractions as low as 


0.1%-1% can be detected [9, 10]. It is from such a detec¬ 
tion that the value of T c can be inferred experimentally. 
This typical range of detection (0.1%-1%) is represented 
by the green area of Fig. 2. We see here that the conden¬ 
sation temperatures calculated at first order and with 
our approach [4] (presented as star symbols in Fig. 2) are 
both located in this windows or in its immediate vicin¬ 
ity. Our generalized formula (solid red line with stars in 
Fig. 2) yields a condensation temperature which defines 
rather well the very beginning of the condensation pro¬ 
cess, when the condensate fraction does not yet vary too 
much with the temperature. It is fully compatible with 
typical experimental determinations of T c . 

Finally, in the last part of the comment, J. Noronha 
discusses different results of the literature obtained with 
power-law traps higher than cubic without settling the 
conflict which exists between these different results in 
these particular types of traps. We note that such a 
debate has also existed for several decades concerning 
the influence of atomic correlations on the condensation 
temperature [7] before this question was finally settled in 
1999 [11]. Answering this question for finite size effects 
requires additional investigations. 
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